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^ . Abstract 

The properties of the ground state wavefunctional in gluodynamics responsible for 
^ ■ confinement are considered. It is shown that confinement arises due to the generation 

of a mass gap in the averaging over the gauge group which is necessary to ensure the 
gauge invariance of the ground state. The string tension is calculated in assumption 
of a particular infrared behavior of the vacuum wavefunctional. 



1 Introduction 



It is generally accepted that the confinement of quarks can be explained as a consequence 
of the linear attractive potential arising between quark and antiquark at large distances due 
to strong vacuum fluctuations of the gluon fields. The origin of these strong fluctuations is 
a complicated structure of the vacuum in QCD. In a literal sense, the vacuum is the ground 
state of the Hamiltonian. In this paper we consider the concrete properties of the vacuum 
wavefunction in gluodynamics which are responsible for confinement. 

In principle, the ground state wavefunctional in QCD can be obtained by solving Schro- 
dinger equation in some reasonable approximation. The crucial point in any practical ap- 
proach to this problem is the maintenance of the gauge invariance. The method of imposing 
gauge invariance condition by integration over all gauges |], appears to be most convenient 
for our purposes. This method is familiar in the gauge theories at finite temperature [[|. [3| . 
It turns out that the averaging over gauges provides a simple way to distinguish between 
Coulomb and confining behavior of the quark-antiquark potential. Confinement arises due 
to the absence of the long-range correlations in the integration over the gauge group. The 
qualitative explanation of the relationship between confinement and the generation of a mass 
gap is given in Ref. J3|. 



2 Ground state wavefunctional and gauge invariance 

The canonical variables in pure gluodynamics (in the temporal gauge A = 0) are gauge 
potentials Af(x) and electric fields E A (x): 

[At(x),E*(y))=i6 AB 6 ij 6(x-y). (2.1) 

We consider gauge group SU (N) and sometimes use matrix notations for gauge potentials: 
Ai = A A T A , where T A are traceless anti-Hermitean generators of SU(N) normalized by 
tr T A T B = -8 AB /2. 

The Hamiltonian of the Yang-Mills theory is 

H = Jd*x (\E A E A + \F A F A y (2.2) 

where F A = d%A A — djAf + gf ABC Af A^ '. Apart from the Schrodinger equation, the ground 
state wavefunction satisfies also the Gauss law constraint: 



DiE A m = 0. (2.3) 

The covariant derivative Di acts in the adjoint representation: D AB = 5 AB di + gf ACB Af. 

In the "coordinate" representation, the electric field operators act as the variational 
derivatives, E A (x) = —i5/5A A (x), and the wavefunction \1/ is a functional of the gauge 
potentials. The operator on the left hand side of Eq. fl2.3p generates gauge transformations 

Ai — ► Af — fi t (Ai + - di] Q, (2.4) 
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so the Gauss law ( |2.3| ) is the infinitesimal form of the gauge invariance condition: 

ty[A n }=^[A]. (2.5) 

Strictly speaking, this is true only for topologically trivial gauge transformations - the wave- 
functional of the (9- vacuum is invariant only up to a phase factor e lke , where k is a winding 
number of the gauge transformation, but for simplicity we assume that 9 = 0, then Eq. fl2.5| ) 
holds for any Q(x). 

In principle, any functional of A can be made gauge invariant by projection on the 
subspace of states obeying the Gauss law. The projection can be realized by averaging 
over the gauge group 0. So, the following ansatz for the ground state wavefunctional 
automatically satisfies ( |2.5| ): 

V[A}= J[DU}e- s[AU \ (2.6) 

The invariance of this expression follows simply from the fact that a gauge transformation 
Ai — > Af can be absorbed by the change of integration variables: U — > VtfU. 

This method of imposing gauge invariance was used in the variational approach of Ref. [|TJ. 
The variational ansatz of Ref. |T[ consists of restricting to the purely quadratic action 

S [A] = \J d 3 xd 3 yA?(x)K(x - y)A?(y). (2.7) 

with some special choice of the coefficient function K(x — y). We consider more general 
situation when the action S[A] can contain arbitrary powers of gauge potentials. However, 
it is clear that the functional S[A] is defined ambiguously. Different choices of S may lead 
to one and the same \I/ and it is convenient to get rid of this arbitrariness demanding the 
quadratic part of S[A] to coincide with So [A]. The color and the Lorentz structure of the 
quadratic form in (|2.7|) is dictated by rotational, translational and global gauge symmetries. 
The general principles allow also a longitudinal term, but it can be always removed without 
change of the wavefunctional, as shown in Appendix. 

The introduction of the "bare" wavefunctional exp(— S[A]), which is not gauge invariant 
and requires the averaging over gauges, nevertheless, has several advantages. There is a 
hope that the bare wavefunctional looks rather simple in some reasonable approximation [|TJ, 
while after the integration over the gauge group one obtains sufficiently complicated gauge 
invariant functional. The possibility to impose the gauge invariance condition exactly on an 
approximate wavefunctional of a simple form makes this method convenient for variational 
calculations JIJ . We shall show that the construction of states containing static color charges 
is straightforward in this framework, which allows to formulate a relatively simple criterion 
of confinement. 

Due to the gauge invariance the normalization integral (^|^) contains a group volume 
factor. This factor is easily extracted in the representation ( |2.6| ): 

(*|tf>= / [DU][DU'][dA] e -S[A u ]-S[AV'] = J [£,£/] [(L4 ] e -SlAV]-S[A) f (2 _ g) 
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where U' is removed by the change of integration variables: — > Af , U — ► U'U. The same 
trick works for the expectation value of any gauge invariant operator PJ. In particular, 



= J [DU][dA] J d 3 x 



'1 PS[A] 1 5S[A) SS[A] 1 " 
2 5AHM 2 5 A? 5 A? 4 4J 4J , 



-S[A U ]-S[A] 



Thus, the vacuum energy is determined by the average of the operator 



R[A] 



d x - 



1 6 2 S[A\ 



2 5 AH At 



1 6S[A] 5S[A] 1 , / 

2 5^ 5^ 4 ij ij . 



(2.9) 
(2.10) 



in the statistical ensemble which is defined by the partition function ( |2.8|) . The same result is 
obtained by averaging of i?^ 17 ] - it arises when the variational derivatives in the Hamiltonian 



act not on the right but on the left. It is natural to preserve the symmetry of (2.8) under 
the interchange of A with Af and to average the symmetric combination iLRLA* 7 ] +i2L4]J. 
Since the vacuum energy is proportional to the volume, it is also natural to add this operator 
to the action and to consider the partition function 



[DU}[dA] e 



-S[A U ]-S[A]- 



-| R[A U ] 



■§R[A] 



(2.11) 



The physical parameters characterizing the QCD vacuum are simply related to the thermo- 
dynamic quantities in the statistical system defined by this partition function. At least, this 
is true for the gluon condensate and the string tension, as shown below. 

The vacuum energy in QCD is equal to the derivative of the free energy In Z with respect 
to A: 



E vac = N ' , ' ' = — \nZ . (2.12) 

It is worth mentioning that the perturbation in ( |2.11| ) softens the short- distance behavior, 
since R is constructed from dimension four operators. So, the parameter A, which has the 
dimension of length, can be regarded also as an UV cutoff. 

The regularized energy density, e vac = E v3iC /V — UV divergent terms, can be expressed 
through the gluon condensate ||: 

/3(a s 



4a, 







pA pAfiv 







(2.13) 



where a s = g 2 and (3(a s ) = — llNa 2 /6n + . . .. The gluon condensate is known to be 
positive. Hence, the regularized vacuum energy density is negative. Therefore, we expect 
that R is negative definite operator after the subtraction of a field independent divergent 
constant ~VK(0) = V J ~ zero point energy which comes from the second variation 

of the quadratic term in S[A]. This explains why the sign before A in Eq. ( |2.11| ) is positive. 



3 Charged states and confinement 



In the presence of charges, the Gauss law (|2.3| ) acquires a right hand side proportional to 
the charge density. We consider the case of two static charged sources placed at the points 
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X\ and x 2 with quark and antiquark quantum numbers. The wavefunctional of such state 
transforms in N (g) N representation of SU(N). The Gauss law constraint for it is 

D t Ef(x) V ab = ig(S(x - x 2 )T b \5 aa , - 5{x - x^T^)*^,, (3.14) 

which implies the following transformation law: 

■V ab [A u ;x 1 ,x 2 \ =rt aa ,{x 1 )n Vh {x 2 )'$ a , v [A-x u x 2 \. (3.15) 

The construction of the state ^^[A; x±, x 2 ) requires operators in the fundamental rep- 
resentation. These operators are necessarily nonlocal in terms of the gauge fields, since 
gluons transform in the adjoint representation. The integral formula ( |2.6|) for the wavefunc- 
tional, in fact, introduces the fields in the fundamental representation, since U transforms as 
U — > WU. This circumstance allows to construct the states with quark quantum numbers 
entirely from local operators at a price of the subsequent integration over all gauges. The 
state obeying (|3.15|) is constructed as follows: 

^ ab [A;x 1} x 2 } = J[DU}U ac ( Xl )U j cb (x 2 )e- s[AU] . (3.16) 

The crucial point here is that the action S in bulk is the same as in the vacuum sector. The 
reasoning is based on the variational arguments. 

The energy of the state ( p,16|) consists of the vacuum energy and quark-antiquark inter- 



action potential: 

#vac + V(X1 ~X 2 ). (3.17) 



The vacuum energy is proportional to the volume, while the interaction potential remains 
finite in the infinite volume limit. Suppose we are trying to minimize the energies of the 
states \1/ and ^ ab . Then 0(V) terms in the variational equations are the same in both cases. 
Since variational equations completely determine the ground state in each sector, the action 
in Eq. ( |3.16| ) coincides with that in Eq. ( |2.6| ) in bulk. In principle, in the charged sector, 
some extra terms localized near the points X\ and x 2 can arise. It is more natural to include 
these terms in the definition of the operators U ac (xi) and Ul b (x 2 ), which create external 
charges. So, strictly speaking, for the state ( |3.16|) to be an eigenfunction of the Hamiltonian, 
the operators U ac (xi) and Ul h (x 2 ) should be replaced by some dressed ones U ac (xi) and 
lA[ b {x 2 ) with the same quantum numbers. The particular form of these dressed operators is 
determined by 0(1) terms in the variational equations. However, in what follows we do not 
distinguish between U and U, since we are interested in a large distance behavior, which is 
determined by the bulk theory and thus is independent of the explicit form of the operators. 

The energy of the state ( |3.16| ) can be expressed in terms of correlation functions in the 
statistical system defined by the partition function (|2.11| ). Analogously to Eqs. ( ^8|) , ( p.9|) 
we find: 

(tfatltfaft) = / [DU\[dA\ trC/(z 1 )C/ t (x 2 )e- 5 ^-^, (3.18) 



^ ah \H\* ab ) = J[DU][dA] tr U{x x )U\x 2 )^(R[A u \ + R[A]) e ~ s ^- s ^. (3.19) 
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These expressions together with Eq. fl2.12p lead to the following representation for the qq 
potential: 

V( Xl -x 2 )= ^-ln(tTU(x 1 )U\x 2 ) 



where the average is defined as usual: 

(0) = ±f[DV)[dA)Oe- 



A=0 



-r[A,u] 



All subsequent consideration is based on this result 



S[A U ]+S[A]-±R[A U ]-±R[A]. 



(3.20) 

(3.21) 
(3.22) 



The averaging over the gauge group in ( 3.21|) may or may not produce a mass gap. 
Consider first the former case. Then the two-point correlator in Eq. ( |3.20 ) falls exponentially 
at large distances: 

tr U{x l )U\x 2 )) = ^ e ~ mr + . . . , (3.23) 



where r = \x\ — x 2 \, m is a mass gap and rj is related to an anomalous dimension of 
the operator U. Note that the constant term in the large distance asymptotics ( |3.23|) is 
forbidden by the symmetry U(x) — > U(x)h, which can not be broken spontaneously, since 
this symmetry follows from the global color invariance of the action S[A}. 

Substituting ( |3.23| ) in Eq. ( |3.20| ) we find that qq potential grows linearly at large distances: 



V{r) 



with 



a 



or + . . . 



dm 



(3.24) 



(3.25) 



A=0 



Thus, the generation of a mass gap in the averaging over gauges leads to the confinement of 
charges in the fundamental representation. It is worth mentioning that the mass gap, and 
thus the string tension a, is determined by the bulk theory. This means that the confinement 
actually is the property of the ground state and not of the operators creating external charged 
sources. 

The screening length in Eq. ( |3.23| ) depends only on the quantum numbers of the operators 
on the left hand side. This justifies the use of the bare operators U instead of the dressed 
ones U in the study of the large distance behavior of the qq potential. The above arguments 
are essentially based on the absence of local operators with quark quantum numbers in 
the pure gluodynamics and they fail for the states with adjoint charges. In this case, the 
state ( |3.16| ) (with the matrices U and TP taken in the adjoint representation) inevitably 
mixes with the ones of the form O a (x\)O b (x 2 )^/ , where O a are some local operators, for 
example, O a = F t A . The energy of such states, roughly speaking, is determined by Eq. ( |3.20| ) 

with (trU(xi)W(x 2 )) replaced by (0 2 (xi)0 2 (x 2 )) . The large distance asymptotics of this 
correlator is determined by an expectation value of O 2 , which is generally nonzero. It is 
easy to see that appearance of the constant term leads to the screening of the interaction 
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potential: V^y(r) ~ e~ mr . However, the mixing with the states composed from purely 
gluonic operators can be small, then the potential for adjoint charges will exhibit a linear 
growth at moderate distances. 

Another remark concerns the anomalous dimension rj entering Eq. (|3.23|) . It gives rise to 
the logarithmic term in the qq potential, unless 



dr] 



0. 



(3.26) 



A=0 



It is generally assumed that subleading terms in the qq potential decrease at infinity, so one 
may expect that Eq. ( |3.26| ) holds for the ground state in gluodynamics. 

Now we turn to the case when the mass gap is absent - m = in Eq. (|3.23|) . Assuming 
the validity of ( |3.26| ) we find that the leading term in the large distance asymptotics of 
the pair correlator is insufficient to determine the behavior of the qq interaction potential. 
Retaining the next term: 



t lU {x 1 )U^x 2 )) = -(l -- + ... 



we obtain the Coulomb law: 



Q 2 



V(r) = const h 



Q 2 



An 



dB 



(3.27) 



(3.28) 



A=0 



Therefore, the presence of long-range correlations in the averaging over gauges leads to the 
Coulomb interaction of static charges. Unfortunately, there is no symmetry that distin- 
guish between Coulomb and confinement phases, in contrast to the finite temperature gauge 
theories g ||. 



4 Examples 

Below the general formalism of the previous section is illustrated by two examples. One 
is £7(1) gauge theory. Another refers directly to gluodynamics. We show that minor as- 
sumptions about an infrared behavior of the wavefunctional allow to calculate the string 
tension. 



4.1 QED 

In the Abelian theory U = e w , Af = A^ + and the action S[A] is quadratic, i.e. 
coincides with ( |2.7| ), and the averaging over gauges reduces to Gaussian integrals. In the 
vacuum sector, the integration over (p in Eq. Q2.6Q replaces A{ by Af- = (5ij — Af 

iff [A] = exp (-- J d 3 xd 3 yA^{x)K{x - y)A^{y)^j . (4.29) 
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This is an exact ground state wavefunctional of the free electromagnetic field, provided that 
K = y/—d 2 (in the momentum space K(p) = \p\). The integration over the gauge group for 
the state fl3.16| ) containing charged sources also can be easily performed yielding: 

V[A;x 1 ,x 2 \ = e i9V(xi) e- i9Vix2) y[A\, (4.30) 

where 

V(x) = J d 3 y -L( x - y)d i A i (y). (4.31) 

The operators e igV ^ were introduced by Dirac ||. The states of type Q4.30] ) created by 
Dirac charge operators behave properly under gauge transformations and are the eigenstates 
of the Hamiltonian of the free electromagnetic field. The wavefunctional ( p.!6|) provides a 
natural nonabelian generalization of this construction. 

Dropping an irrelevant constant term we get for the action (|3.22|) in the Abelian theory: 

T = \J d 3 xd 3 yA t (x) (k(x - y) + ^K 2 (x-y) ) j A(y) + (A t <— At+fycp)-^ J d 3 xF tJ F tJ . 

(4.32) 

Here K 2 (x—y) denotes the kernel of the operator K 2 . The gauge potentials can be integrated 
out by the change of variables: 

Ai = A\ - i d i( p. (4.33) 
The shifted gauge potentials A\ decouple and we are left with the effective action for <p: 

r cfr = \ J d 3 xd 3 yd lV (x) (k(x -y) + \ K 2 (x - y)\ d iV {y). (4.34) 

So, for the two-point correlator we get: 

' e igtp{xi) e -igtp{x 2 )\j _ e g 2 {D(x 1 -x 2 )~D(0)) ^ (4.35) 

where D is the propagator of the field (p: 

d " = z t{ k+ - 2 k2 )- ( 4 - 36 ) 

Differentiating the propagator with respect to A we find: 

1 



3D 



A=0 U 



(4.37) 



Hence, according to Eq. (|3.20|) , we obtain the Coulomb law for the interaction potential: 

V(xi - x 2 ) = -g 2 - ^2) - = + self-energy. (4.38) 
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It is interesting that the explicit form of the coefficient function K was not used anywhere in 



the derivation of this result. Nevertheless, it is useful to substitute K = -J — d 2 in Eq. Q4.36D 
in order to compare ( 4.35 ) with ( |3.27| ). After some calculations we obtain: 



2r ■ 2r „• 2r ■ 2r N 



2r 
2r 




(4.39) 



where 7 is the Euler constant. Comparing this expression with (|3.27| ) we find that the 
constant C diverges in the limit A — > 0. This divergence reflects the fact that the self-energy 
of the Coulomb charge is infinite, but this infinity is regularized by A. 



4.2 QCD 

Due to the asymptotic freedom, the short distance behavior of the qq potential is determined 
by the perturbation theory. At weak coupling, one can expand U = e 9 * T = 1 + g<& A T A + 
. . .. Only the quadratic term in the action S[A] contribute to the leading order in g and 
the calculations are literally the same as in the Abelian theory. Since the explicit form of 
the coefficient function K is inessential for these calculations, at short distances we get the 
Coulomb potential Q4.38 ) with g 2 multiplied by a group factor (iV 2 — 1)/N. 



It is naturally to assume that at large distances the long wavelength approximation is 



valid and the action ( 3.22 ) can be approximated by its derivative expansion. The structure 
of the derivative expansion for T depends on the behavior of the coefficient function K{jp) at 
small momenta. Again, only the quadratic term in the action is essential, since the remaining 
ones correspond to operators of higher dimension. Here we assume that K(p = 0) 7^ 0, so 

K(p) = M + 0(p 2 ). (4.40) 

Neglecting 0(p 2 ) terms we obtain: 

T = ^Jd 3 x(A?A? + AY A A? A ), (4.41) 

where 

AM 2 

pi = M+—. (4.42) 

The magnetic term (F-) 2 is also irrelevant in this approximation, since it contains extra 
derivatives in comparison to (Af ) 2 . 

The case of purely quadratic action S with the coefficient function K approximated by a 
constant at small momenta was considered in Ref . U . Following this consideration, which is 
based on methods typical for nonlinear sigma-models, we express the string tension in terms 
of M and the gauge coupling constant. 

First, we get rid of the integration over Ai by the change of variables: 

Ai = Ai -1-diUlfl, (4.43) 
2# 
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which is a counterpart of ( J4.33|) . After this, A4 decouples and we get the effective action for 
U: 

T cS = J d 3 x tr (y^diU) 2 = ^- 2 Jd?x todiV%U. (4.44) 

It is convenient to integrate over U(N) instead of SU(N). This can be achieved by adding 
the following term to the action: 

Ti = J d 3 x (trrfdtU) 2 . (4.45) 

Since T c g + I\ is invariant under phase transformations U — > e tip U, it depends only on 
SU(N) variables, for which it reduces to r eff . The last step consists in the replacement of 
the group integral by an integral over complex N x N matrices. The unitarity condition is 
imposed by means of a Lagrange multiplier: 



JL 
2g 2 



J d 3 x tra (U ] U - l) , (4.46) 



where 

a ab = m 2 5 ab + v ab , (4.47) 

and v is anti-Hermitean matrix. 



The crucial point is that (|4.46[) introduces a mass term for U, so the pair correlator falls 
exponentially at large distances as in Eq. ( |3.23| ), which leads to confinement. The mass gap 
m is determined by the equation of motion for a. This equation follows from the effective 
action obtained by integrating out U(x). The gap equation at the same time is equivalent 
to the unitarity condition (U(x)W(x)\ = 1. In the one-loop approximation, we have 



d3p 1 - " (4.48) 



(2tt) 3 p 2 + m 2 2g 2 N 

The left hand side is UV divergent and, in principle, it is necessary to introduce a cutoff in 
order to find the mass gap 0. However, the string tension appears to be cutoff independent, 
since for the derivative of m we obtain: 

dm 2 r d 3 p 1 M 2 

1 (4.49) 



OX J (2tt) 3 (p 2 + m 2 ) 2 Ag 2 N' 



and, according to Eq. (|3.25|) , the string tension is 



dm nM 2 M 2 

° = ~ -dx = ?iv = wv" (4 - 50) 

This result is not universal, it is valid under the condition that K(p = 0) 7^ 0. This 
assumption is compatible with the variational estimates of Ref. [0. However, the arguments 
in favor of K(p = 0) = also exist, and there are some indications @ that the expansion of 
K(p) at small momenta begins with the term of order p 2 . If such behavior actually holds, 
the present treatment of confinement properties should be modified. 



10 



5 Discussion 



Our main result is that, in the Hamiltonian picture, the confinement arises as a consequence 
of the generation of a mass gap in the averaging over gauges. Although it is possible to derive 
the linear potential and even to calculate the string tension under some simple and not very 
restrictive assumptions about the infrared properties of the ground state in gluodynamics, 
the quantitative treatment requires a more detailed knowledge of the vacuum wavefunctional. 
Possibly, the variational techniques |L] can provide a reasonable approximation for the ground 
state in QCD incorporating both asymptotic freedom || and confinement. 

A more detailed quantitative treatment of the ground state in gluodynamics may provide 
insights in some other nonperturbative properties of the low-energy QCD. In particular, the 
usual picture of the qq interaction is associated with a string stretched between quark and 
antiquark. Perhaps, this picture would correspond to a sum-over-path representation for the 
correlator entering Eq. ( |3.20 ): 



t-X (1)=X2 

(trU(x 1 )rf(x 2 ))= [dX] e"*™. (5.51) 

\ ' JX(0)=x 1 



The large distance asymptotics (|3.23|) is naturally interpreted as a saddle point approxima- 
tion for the path integral (|5.51 ) and the functional 7i[X(£)] has a meaning of an effective 
string Hamiltonian. 
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Appendix A Removing arbitrariness in the action 



To demonstrate that the representation (|2.6|) is ambiguous let us consider a set of functionals 
S[A; £], whose dependence on the parameter £ is determined by the equation 



^ = / d*xd 3 yA?(x)W(x - y)diDt B j§^ y (A.l) 



if 

Here W(x — y) is an arbitrary function. Then the wavefunctional 



= J[DU] e - s[AU * ] (A.2) 

is independent of £ up to a normalization factor, which is inessential. The proof is based on 
Schwinger-Dyson equations. 

The measure of integration in ( [A.2 ) is invariant under left shifts: 



5JJ = Uu, 5 W U ] = -ujU\ (A.3) 
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where uj is traceless anti-Hermitean matrix. The twisted gauge potentials Af transform 
under left shifts as 

6 U A? = [A?, u] + ldiV = - D^u. (A.4) 
The Schwinger- Dyson equation follows from the identity 



-S[A l 



= J[DU]5 u (dtA? A (x)e- s ^ 

= lj[DU\ (d i Df AG u J c (x)+d i Af A (x) J d 3 y^(y)Df CB ^^j e 

This equality holds for any uj c (y). Taking uj c {y) = 8 CA W(x — y) and summing over a we 
find: 

(N 2 -l)d 2 W(0)V[A}+ f [DU] f d 3 y cUf A (x) W (x - y) Df AB /fffl e~ s ^ = 0. (A.5) 

J J ' oAj \y) 



Combining this equation with (|A.l ) we obtain: 



^ = -J[DV\ J dVyAf A (x)W(x-y)d l Df AB I J ] ^- ) e"^ = c% (A.6) 

where c= (N 2 - l)Vd 2 W(0). Thus, 

*L4;£] = e^[A;0}. (A.7) 

The above arguments show that the solution of Eq. ( |A.1| ) for any W and at arbitrary 
value of the parameter £ gives essentially the same wavefunctional. The action S can be 
expanded in the power series: 

S[A] = \J MyAt^Kyix - y)Af(y) + .... (A.8) 



Then the equality ( A.l ) is rewritten as an infinite set of equations for coefficient functions. 



The first equation in this set contains only the two-point function and reads: 

dK 

= WdAK kj . (A.9) 

Here W denotes an integral operator with the kernel W(x — y). We can start with the 
general expression compatible with translational and rotational invariance taken as an initial 
condition for Eq. ( |A.9| ): 



K,M=[^-^§ L )K+^K', (A.10) 
where K and K' are the functions of d 2 only. Solving Eq. (|A.9|) we find: 



K tl (() = L - d ~0) K + °& e^K>. (A.ll) 
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Taking 

£W = (d 2 )- l (\nK -\nK') (A.12) 

we get the diagonal coefficient function Kij = SijK, thus proving the assertion used in Sec. ||. 
The extreme case of £W — > +00 corresponds to purely transversal coefficient function, which 
arises if the functional S[A] is gauge invariant. 
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